The Hellberg-Mele Jastrow factor as a variational wave function for the one 

dimensional XXZ model 
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We find the Jastrow factor introduced by Hellberg and Mele in their study of the one dimensional 
t — J model provides an exceedingly good variational description of the one dimensional XXZ model. 
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The one dimensional systems are genuinely strongly 
correlated. This is most clearly indicated by the power 
law behavior of various correlation functions, which are 
termed Luttinger liquid behavior in general Ij . In a Lut- 
tinger liquid system, the powers of the correlation func- 
tions, or, the critical exponents of the system, are con- 
tinues functions of the model parameter. 

The one dimensional XXZ model is typical Luttinger 
liquid system. In this short note, we propose a variational 
description this well known model. The Hamiltonian of 
the model reads 

H = - J2(S^S^+, + Srs^, + ASf S^i), (1) 

i 

in which A is the parameter of anisotropy This model 
is exact soluble in terms of the Bethe Ansatz 2]. For 
A < — 1, the system is in the Ising regime in which the 
ground state is antiferromagnetically ordered. For A > 1, 
the system phase separates into full polarized regions. 
For — 1 < A < 1, which is the most interesting case, 
the system exhibits critical behavior with continuously 
varying critical exponents. For example, the asymptotic 
behavior of the transverse and the longitudinal spin cor- 
relation function are given by 
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in which the critical exponent rj is given by 

1 1 • -1 A 
ri = sm A. 
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The last equation for the critical exponent rj is obtained 
by combining the result of Bethe Ansatz solution and the 
effective field theory based on Abelian Bosonization|3. 
Either the Bethe Ansatz solution or the effective field 
theory alone is not powerful enough to predict the critical 



exponent. In the Bethe Ansatz approach, the wave func- 
tion is so complicated that a direct calculation of the cor- 
relation function is impossible. On the other hand, while 
the essence of the critical correlation is well captured by 
the effective field theory approach, the parameters in the 
theory must be set by hand. 

A good variational wave function should capture simul- 
taneously the short range and the long range correlation 
of the system studied. However, for local Hamiltonian, it 
is difficult to get the correct long range behavior of the 
system by optimizing the energy of a variational wave 
function. Wave functions very close in energy can have 
drastically different long range behavior, provided that 
the system is in the critical regime. 

We find a wave function originally proposed for the one 
dimensional t — J model provides an excellent description 
of both the short rang(energy) and the long range correla- 
tion of the one dimensional XXZ model for —1 < A < 1. 
The wave function is given by 
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in which Xi denotes the coordinates of the up spins in 
the periodic chain of N sites. This form is first proposed 
by Hellberg and Mele to describe the Luttinger liquid 
behavior of the one dimensional t— J model 4] . In terms 
of the t — J model, this function appears as a Jastrow 
factor for the residual charge correlation in front of the 
well known Gutzwiller projected Fermi sea wave function. 
For that purpose, one should reinterpret ^msi{{xi}) as 
a wave function for hard core Boson, in which Xi then 
denotes the coordinates of the charges(or holes) in the 
t — J model. 

The function introduced by Hellberg and Mele is 
also known as the exact ground state of the Suther- 
land model with inverse square interaction in one spatial 
dimension_5j. The Sutherland model reads, 
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where the last term is a generahzed inverse square po- 
tential on a periodic chain of length L. It is found that 
when 



= (v/25TT+l)/2, 



(7) 



^Hivids^i}) is the exact ground state of the Sutherland 
model. 

In our recent work on the variational study of the one 
dimensional t — J model, we find the residual charge 
correlation beyond the Gutzwiller projected wave func- 
tion(GWF) of the model should be described by a XXZ- 
type effective IIamiltonianjy|. Combining this analysis 
and wave function proposed by Hellberg and Mele, one 
quickly realized that the Jastrow factor introduced by 
them should also be a good description of the ground 
state of the one dimensional XXZ model itself. 

The reason for the excellentness of the Hellberg-Mele 
wave function for the one dimensional XXZ model can be 
more directly seen as follows. At A = 1, the XXZ model 
reduces to that of the isotropic spin chain with ferromag- 
netic exchange, whose ground state is the fully polarized 
state. The wave function for the fully polarized state is 
a constant in the Ising basis, which is given exactly by 
^HM({a;i}) with u = Q. At A = 0, the model reduces 
to the one dimensional XX model, which is equivalent 
to the free spinless Fermion model through the Jordan- 
Wigner transformation. In this case, the ground state 
wave function is given exactly by ^HM({a;i}) with v — 1, 
a Slater determinant for the spinless Fermion. Another 
special case is when A = — 1. In this case, the model 
reduces to that of the isotropic spin chain with antifer- 
romagnetic exchange. Although ^Hivida^i}) is no longer 
an exact ground state for this model, it is well known 
that the '^BM{{xi\) with u = 2, namely a Gutzwiller 
projected half filled Fermi sea wave function, provides 
a exceeding good variational description for the model 
consideredf^. For example, the energy calculated from 
^HM({a;i}) is about -0.6921 per bond, which very close to 
the exact value, — In 2. At the same time, the long range 
correlation of the model is also correctly captured by this 
wave function, apart from a logarithmic correction which 
is absent for A > —1. 

The ground state energy and other ground state corre- 
lation of '^BM{{xi}) can be easily calculated by the Vari- 
ational Monte Carlo method. Figure 1 shows the ground 
state energy calculated from ^HM({a:i})- We find the 
relative error in the ground state energy is less than 0.24 
percent ( which is reached at A = — 1) for all values of A. 
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FIG. 1: Variational ground state energy calculated from 
'^HM{{xi\). The solid line shows the exact result obtained 
from Bethe Ansatz solution. The optimized value of the A is 
shown in the inset. 



the variational result agree qualitatively with the exact 
one. Note for A = 1, although the spin rotational sym- 
metric GWF give the correct value of 1 for the critical 
exponent ?/, one find that in the variational description 
in terms of ^mAi^i) the system chooses to break such 
symmetry slightly. We think this failure should be at- 
tributed to the logarithmic coorection at the symmetric 
point A = — 1. 
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The Hellberg- Mele- type wave function can not only 
give precise estimate for the ground state energy, but 
also give qualitatively correct behavior of the critical ex- 
ponent of the model. In [Sj , it is found that the critical 
exponent 77 for ^nM{xi) is simply given by |-. In Figure 2, 
we compare the exact result for the critical exponent with 
that given by the Hellberg-Mele wave function. One find 



FIG. 2: Correlation exponent r) as defined in Eq.(l) as deter- 
mined from the Hellberg-Mele wave function. The solid line 
denotes the exact result ri = 1 — — cos~^ A. 
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Finally, we note that although the Hellberg-Mele wave 



function is good approximation for the S^ 



sector 



of the one dimensional XXZ model, it fails to describe 
the physics of the Sjo^ ^ sector of the same model. In 
p , this is found to be responsible for the failure of the 
Hellberg-Mele-type wave function to describe the Luther- 
Emery phase at small electron density and large J jt of 
the one dimensional t — J model. 
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